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Abstract — A new class of spherical codes is constructed by 
selecting a finite subset of flat tori from a foliation of the unit 



sphere S 



and designing a structured codebook on 



each torus layer. The resulting spherical code can be the image 
of a lattice restricted to a specific hyperbox in R L in each layer. 
Group structure and homogeneity, useful for efficient storage and 
decoding, are inherited from the underlying lattice codebook. A 
systematic method for constructing such codes are presented and, 
as an example, the Leech lattice is used to construct a spherical 
code in R . Upper and lower bounds on the performance, 
the asymptotic packing density and a method for decoding are 
derived. 

Index Terms — Spherical codes, group codes, flat torus, lattices, 
Gaussian channel. 



I. Introduction 

The problem of placing points on the unit Euclidean sphere 
of a given dimension has attracted the attention of engineers, 
mathematicians and scientists and has relevance to many di- 
verse fields of science and engineering. In communication the- 
ory, point sets on the unit sphere are useful for communicating 
over a Gaussian channel and are a natural generalization of 
phase shift keyed signal sets (PSK) to dimensions greater than 
two. The point placement problem in this case is formulated as 
a packing problem in which the objective is to pack as many 
spherical caps of given radius as possible on the sphere. The 
dual to the packing problem is the covering problem, useful 
in facility location, in which the problem is to minimize the 
number of spherical caps of a given radius so that every point 
on the sphere is covered. In quantization, point sets on the 
sphere form a key component in shape-gain vector quantizers 
0. Point sets on the sphere with special properties are known 
as spherical i-designs. An overview on spherical codes and 
its properties can be found in lfl4ll . and two families of 
asymptotically dense codes are presented in [8|, |9|. Lists of 
good spherical packings, coverings and designs can be found 
online at lfL3l 

In this paper, we describe a new method for constructing 
spherical codes for the communication problem. While it is 
important to maximize the packing density, additional practical 
considerations such as storage and easy decoding are also 
important. The codes presented here have low construction 
and decoding complexity and, for not asymptotically small 

C. Torezzan is with the School of Applied Sciences, University of 
Campinas, Brazil (e-mail: cristiano.torezzan@fca.unicamp.br. S.I.R. Costa 
is with the Institute of Mathematics, University of Campinas, Brazil (e- 
mail: sueli@ime.unicamp.br. V. A. Vaimshampayan is with AT&T Shannon 
Laboratory, Florham Park, NJ, USA (e-mail:vinay@research. att.com). 

This work was supported in part by Sao Paulo Research Foundation 
FAPESP and National Counsel of Technological and Scientific Development 
- CNPq - Brazil. 

The material in this paper was presented in part at the IEEE International 
Symposium on Information Theory, Seoul, Korea, 2009. 



distances, have comparable performance to some well known 
apple peeling J6), wrapped [8| and laminated [9] codes. 

The paper is organized as follows. The flat tori foliation 
of the sphere is introduced in Sec. [TT] Sec. Ill describes our 



proposal of Torus Layes Spherical Codes (TLSC) which is 



IV 



based on a foliation of the unit sphere in R 2L . In Sec 
an example of TLSC in i? 4 , which is cyclic in each layer is 
presented with comparisons to some well known constructions. 
In Sec. |V]it is described how lattices with good packing density 
in R L can be used to construct TLSC in R 2L . As an example 
the Leech Lattice is used to construct a spherical code in R 48 . 
Upper and lower bounds on the performance and asymptotic 
packing density of the TLSC are derived Sec. VI Finally a 



decoding method for our codes is described in Sec. VII 



II. Foliation of the Sphere by Flat Tori 

The unit sphere S' 2i " 1 C M 2L can be foliated with flat tori 
(Clifford Tori) (T), as follows. For each unit radius vector 
c = (ci, c 2 , .., c L ) £ R L , Ci > 0, and u — (ui,u 2 , ■ ■ ■ , u L ) G 
R L , let $ c : R L R 2L be defined by 



■r / \ I I Hi . Ul 

<P c (u) = ci cos — , sin — 

Ci Ci 



Ul . UL 

, cl cos — , sin — 
cl c L 



(1) 

The image of $ c is the torus T c , a flat i-dimensional surface 
on the unit sphere S 2L ~ 1 . Q c is an embedding of the flat torus 
T c , generated by the hyperbox: 



V c = {u& R L ;0 <u % < 2-KCi}, 1 < i < L. 



(2) 



Note also that each vector of S 2L 1 belongs to one, and 
only one, of these flat tori, some of which may be degenerat^] 




Fig. 1. Ilustration of a torus layer spherical code in dimension four, projected 
in M 3 . The codewords belong to the surface of flat tori. 

'A degenerate torus T c is one for which some c, = 0, i = 1, 2, . . . , L, 
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The Gaussian curvature of a torus T c is zero [1 J and T c can 
be cut and flattened into an L-dimensional box, V c , just as a 
cylinder in M 3 can be cut and flattened into a 2-dimensional 
rectangle. Since the inner product (d$> c /dui, d$ c /duj) — Sij, 
the application $ c is a local isometry, which means that any 
measure of length, area and i-dimensional volume on T e is 
the same of the corresponding pre-image in the L-dimensional 
hyperbox V c 

We say that the family of flat tori T c and their degenerations, 
with c = (cx,C2, ■■,Cl), \\c\\ — 1, Ci > 0, defined above is a 



2L-1 



C 



I>2L 



foliation on the unit sphere of S 

Let Tf, and T c be two flat tori, defined by unit vectors b and 
c with non vanishing coordinates. We can assert that: 

Proposition 1: The minimum distance between two points 
on these flat tori is 

\ 1/2 



d(T c ,T b ) 




Proof: This follows from observing that 



D 



\^ C (U) - <P C (V)\\ 2 > \\c-b\\ 



(3) 



(4) 



with equality holding if, and only if. 

V4 _ U 

a t 



Vi 



Note that Proposition [JJ can also be extended to degenerate 
tori by replacing c, ^cos I 1 , sin by (0, 0) if c,- = injlj 
The distance between two points on the same torus T c given 

by 



d($ e (u),$ c (v)) = 2JJ2 c h™ 2 ( 



2c, 



) 



is bounded in terms of \\u — v\\ by the following proposition. 

Proposition 2: Let c= (ci, C2, .., Cj,), ||c|| = 1, and let 
u,v e V c - Let A = \\u-v\\ and 5 = || $<•(«) - *c(«)||- Then 

„;„ A .A 



-A < 



A 

2^ 



2c<: A < (5 < ^A < A 



A 

2 



where ^ = argmin(cj). 



Proof: This proposition is an extension of a re- 
sult presented in lfT6l and its proof follows similar 
arguments. We can show that for fixed A, the mi- 
nimum and the maximum distortion, which correspond 
to maximum and minimum bending, occur whenever 
u — v = Ac and u — v = Ae^ respectively, where ej denotes 
the canonical unit basis vector which is nonzero only along 
the jth coordinate. ■ 



III. Torus Layer Spherical Codes in 



t>2L 



Our goal is to construct a spherical code in E with 
minimum distance equal to a given value d. We denote such 
a code by TLSC{2L,d). 

Before presenting a formal construction technique, we de- 
scribe the main idea. Given a distance d g (0, V2], we firstly 
define a finite set of tori on S 2L ~ 1 such that the minimum 




Vc 



2nc 




Fig. 2. Ilustration of the construction of a four dimensional torus layer 
spherical code. 



distance, according to Proposition [TJ between any two of these 
tori is greater than d. Then, for each one of these tori, a finite 
set of points is chosen in M. L such that the distance between 
any two points, when embedded in M 2i by the standard 
parametrization ([TJi, is greater than d. This set of points may 
belong to a L— dimensional lattice, restricted to the hyperbox 
^ or to any other suitably chosen set. The TLSC(2L,d) is 
the union of the images under ([TJ of each finite sets of points, 
one for each torus. Figure [2] illustrates the construction of a 
TLSC{A,d). 

Note that the pre-image of the points in a single layer of 
TLSC(2L, d) lie inside an L— dimensional box and hence we 
are working in half of the code dimension. For not that small 
values of d, our codes compare favorably in terms of code 
size with previous codes (5), (8) and J5] (see Tables IV and 
[TTT^. In addition, the group structure of our code in each layer 
allows efficient storage and decoding. 
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A. The construction of Torus Layer Spherical Codes 

Let L > 2 and d e (0,\/2]. Let SC{L,d) be an 
L— dimensional spherical code with minimum distance greater 
than d. The code TLSC(2L, d) is constructed in two steps as 
follows: 

(i) Select the points in the SC(L, d) which have only 
nonnegative coordinates. This sub-code is denoted by 

SC(L, d)+ = {ce SC(L, d), c t > 0, 1 < i < L} . 

Each point c £ SC(L, d) + defines a fiat torus T c in the 
unit sphere in M. 2L and hence a hyperbox V c , according 
to (|2J.. 

(ii) For each torus T c , defined by SC(L, d) + , look for the 
largest set of points Yr c C V c such that 

||$c(j/) - $e(*)|| >d Vx,yeY Tc . 

The performance of a TLSC is directly related to the 
methods used for constructing SC(L,d) + and Yt c - In (i) 
we must choose a SC(L,d) + with good density and, if 
possible, some symmetries. For this purpose we may use any 
suitable spherical code, e.g., a L-dimensional TLSC, or some 
other known structured spherical codes, such as wrapped [8|, 
laminated [9| or apple peeling |6| spherical codes. We could 
also use a non structured spherical code, e.g. one of the codes 
listed at [ 1 3 1 . Since the cardinality of the set SC(L,d)+ is 
much smaller than the final code, unstructured spherical codes 
are also acceptable. For (ii) a good option is to consider lattice 
points inside each hyperbox V c . Through the maps $ c they 
generate group codes in each torus layer ATI . In the next 
Sections we present examples of this construction. 

A two step-TLSC(2L + l,d) can also be constructed in 
odd dimensions by first slicing the unit sphere S 2L C K 2L+1 
through hyperplanes perpendicular to the canonical vector 
S2L+1, such that the minimum distance between two 
hyperplanes is at least d. Table IV shows a comparison 



where 



between TLSC (5, d) and the apple peeling codes presented 
in 1 6 1 for the same distance. 



IV. A PIECEWISE CYCLIC FOUR DIMENSIONAL TLSC 

In order to clarify the technique and present an example, we 
construct a torus layer spherical code TLSC(A,d). For each 
layer we will design a cyclic group code so that the resulting 
code will be a piecewise cyclic spherical code. 

Step (i) in this construction is to choose a good spherical 
code in M 2 . Since the best spherical code in M 2 with minimum 
distance d is unique (up to rotation) and is symmetric (the code 
is the set of vertices of a regular polygon inscribed in the unit 
circle), the only design choice is to determine a good rotation 
for points in the positive quadrant. 

Our approach is to select points located symmetrically in 
relation to the line with unit slope. Thus 

SC(2,d)+ = {(cos(a ±j ),sin(a ±J )),0 < a ±j < |}, (5) 



7T .id 

a±j = — ± (2j — 1) arcsm 



4 

1 < j < 



tt—2 arcsin (d/2) 
8 arcsin (d/2) 



(6) 



«2 = J + 3 arcsin d/2 

cii = J + arcsin d/2 




Fig. 3. SC(2, symmetric in relation to y = x 

Note that each flattened torus T a± . is a rectangle with sides 
of length 2ircos(a±j) and 2-7T sir^a-tj). 

In step (ii), we fill each torus T a±j searching for the largest 
cyclic group code, with minimum distance d and initial vector 
x 0±J = (cos(a±j),0,sm(a±j),0). 

A cyclic group code of order M in R A is an orbit of an initial 
unit vector under action of the cyclic group of orthogonal 
matrices generated by a (2 x 2)-block diagonal matrix 

Gj = 
I cos( 



M 



\ 



sinC 2 ^) 

M I 






cos(- 





cos(^) 
-sin(^) 






cos(^ 



■) 

-) / 



where gcd(gj 1 , gj 2 ) — 1. The pre-image of this cyclic group 
code by $ c is a lattice in R 2 ITT1 . 

The search for the largest group code for each torus T a ., 
attained by an L-dimensional lattice can be accomplished 
based on [15] and a simplified algorithm is described here. 

This problem is, loosely, a dual to the initial vector problem 
(IVP) fl2l . 10, which is a classic problem in group codes. In 
IVP is given a group and the problem is to find a unit vector in 
order to maximize the minimum distance in its orbit. Here we 
have the initial vector xq, (which defines the torus T a± .) and 



wish to find the largest cyclic group such that the minimum 
distance in its orbit is, at least, a previously fixed value d. 

A. An example: TLSC{4,0.3) 

We now illustrate the construction of a quasi-cyclic torus 
layer spherical code for n = 4 and d = 0.3. 
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Input: d, a 

Output: Generators: {gj 1 ,gj 2 } 
xq = (cos(a), 0, sin(a), 0); 
7r 2 cos(a) sin(a) 

L2 % /3arcsin(|) 2 _ 
continue = 1; 
while continue do 
for g n = 1 to [f 
for <? J2 = 1 to 



do 

fj do 
if gcd(g h ,g h ) = 1 then 

<i ~ min (G,) l a;o — h 

_i<»<LfJ 

if d > d then 

Return {gj^gj.,}; 
continue — 0; 
Stop; 

end 
end 



end 

M -- 



end 

M 



1; 



end 

Algorithm 1: Algorithm to search for the best cyclic group 



code in R 4 for a given initial vector x 0j . 

• From [6] we get a + i = 0.935966, a+2 = 1-2371, and 
a +3 — 1.53824, which define the points of SC(L,d) + 
above the line y = x, according (|5j. 

• For each torus T a+ , we have found the largest cyclic 
group code using the algorithm [T] The result is: 



a 


cos(a) 


sin(«) 


dmin 


M 


(Sa,3i2) 


0.935966 


0.593041 


0.805173 


0.30225 


233 


{1,98} 


1.237103 


0.327535 


0.944839 


0.301406 


146 


{22,1} 


1.538240 


0.032551 


0.99947 


0.312869 


20 


{0,1} 



TABLE I 



Part 1 of TLSC(i, 0.3): tori above the slope line. 



• Finally, for each torus T a+j we consider the symmetric 
layer T a _ ., just interchanging the coordinates. 





a 


sin(a) 


cos(q) 


dmin 


M 


Generator 





634829 


0.805173 


0.593041 


0.30225 


233 


{98,1} 





333694 


0.944839 


0.327535 


0.301406 


146 


{1,22} 





032559 


0.99947 


0.032551 


0.312869 


20 


{1.0} 



TABLE II 



Part 2 of TLSC(i, 0.3): tori below the slope line. 



The resulting code TLSC(4,0.3) has 6 layers, pairwise 
symmetric, with 20, 146, 233, 233, 146, 20 points respectively 
and thus the entire code has 798 points. 

Due to the symmetry and group structure of this code, in 
order to store all 798 codewords in this TLSC(A, 0.3) is only 
required columns 1, 5 and 6 of table [I] The constructiveness of 
the codewords is a good aspect of the Torus Layers Spherical 
Codes when compared with several other known construction 
of spherical codes. 



In the Table |TTTJ we compare torus layer spherical codes to 
three other known spherical codes: apple-peeling [6], wrapped 
JS) and laminated (9J, at various minimum distances d. 



d 


TLSC(4,d) 


apple-peeling 


wrapped 


laminated 


0.5 


172 


136 






0.4 


308 


268 


* 


* 


0.3 


798 


676 


* 


* 


0.2 


2,718 


2,348 




* 


0.1 


22,406 


19,364 


17,198 


16,976 


0.01 


2.27 xl0 Y 


1.97 xl0 Y 


2.31 xl0 Y 


2.31 xl0 Y 



TABLE III 

Four-dimensional code sizes at various minimum distances. (*): 
unknown values. 



Using four dimensional codes and successive slices of the 
S 4 C M 5 by hiperplanes we constructed some TLSC(5,d). Al- 
though S 4 is not foliated by flat tori, the codes constructed in 
layers of tori in S 3 and lifted to S 4 exhibit good performance 
as illustred in Table [TV] 



d 


TLSC(5,d) 


APC(5,d) 


0.8 


48 


48 


0.7 


98 


64 


0.6 


196 


160 


0.5 


374 


336 


0.4 


872 


872 


0.3 


3,232 


2,960 


0.2 


17,140 


15,424 


0.1 


296,426 


256,760 


0.05 


4,824,018 


4,164,152 



TABLE IV 

Comparison between 5-dimensional torus layer and apple 
peeling spherical codes at various minimum distances 



V. Orthogonal sublattices and piecewise 

COMMUTATIVE GROUP CODES 

In this section we discuss how dense lattices in R L with 
good packing density can be used to construct a well structured 
TLSC in R 2L . This is done by considering orthogonal L- 
dimensional sublattices in step (ii) of our construction. 

We first complete the step (i) of our construction for a given 
minimum distance d where each point c G SC(d, L) + C 
5 L_1 defines precisely the lengths of an orthogonal hyperbox 
V c C K L . For a chosen L-dimensional lattice A, we look for 
an orthogonal sublattice Ai such that the fundamental region 
F of Ai approaches V c . Then Ai should be scaled to a lattice 
Ai in order to get V c as its fundamental region. 

The lattice points of Ai inside V c are identified with the 
quotient A/Ai and inherit the associated group structure, since 
A/Ai is isomorphic to A/Ai. The image by Q c of these points 
defines a commutative group code in M 2L , with initial vector 
defined by vector c € SC(d, L) + |Qj], iTBl . 

In addition, let B and B\ be the generator matrices of the 
lattices A and Ai respectively and let Q be an integer matrix 
such that Bi = BQ. Then the characterization and the set 
of generators of the group A/Ai can be obtained from the 
standard Smith Normal Form of Q 0, H3. 
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This strategy can be recursively applied to all flat torus 
layers defined by SC(d, L) + to get a piecewise commutative 
TLSC in R 2L , i.e. a spherical code constituted of layers of 
commutative group codes. 

A. A 48-dimensional TLSC from the Leech Lattice 

To illustrate the construction described above, we present 
next a piecewise commutative group code in i? 48 with mini- 
mum distance d = 0.1, designed from an orthogonal sublattice 
of the Leech Lattice. Although the number of points in this 
code is of order 10 34 the construction is quite simple and 
do not need the storage of the points. In addition, there is 
a fashion labeling for all points in this code induced by the 
commutative group code in each layer. 

In order to simplify the step (i) and focus our attention in 
step (ii), we start from a set of 24 points in the unit sphere 
S 23 consisting of all permutations of the vector 

V ' v / 23+^ 

where t(d) > is chosen in order to guarantee the desired 
minimum distance d between any two permutation. 

In this example, to design a code with minimum distance 
d = 0.1, we get t(0.1) = 1.35234 and 



c(1.35234) °£ (0.271399, 0.200688, • • • , 0.200688) G S 23 

Each permutation of c(t) defines a flat torus on the surface 
of S C M 48 which can be flattened into a 24-dimensional 
hyperbox V c - Therefore, the step (i) of our construction is 
done. 

The next step is to find a discrete set of points inside each 
hiperbox defined above. Since the 24 hiperbox differs only by 
rotations (or interchange of coordinates) it is enough to solve 
this problem for one of those tori and we show next how to 
use the Leech lattice to solve that. 

Consider a standard Leech Lattice rescaled by a factor 
/3(d) — 0.10187 to assure minimum distance d = 0.1 in i? 48 , 
after apply the function 3> c (t) (according Proposition 12b. 

Let B be the generator matrix of the scaled Leecn lattice 
A.2Ap with minimum distance j3. Since the Leech lattice 
contains the sublattice 4Z 24 , we can determine a factor a such 
that the lattice generated by al, where I is the identity matrix 
of order 24, is an orthogonal sublattice of A 2 4^. 

Now, let v(t) be an integer vector defined by 

2nc(t) l - (3 

Vi = 

a 

Since the length of each edge of the flattened hiperbox V c 
is given by 2irCi(t), each coordinate i of vector v determines 
the maximal number of times vector aei can be placed in 
each canonical direction of V c , i.e. vector v defines the largest 
orthogonal sublattice of A24/3 in V c . In our example 

«(1.35234) = (11,8,8, ••• ,8). 

So the matrix B\ = adiag(w) is a generator matrix of a 
sublattice A 2 4^ C A24/3 which fundamental region approaches 



hiperbox V c - Therefore the number M,;(i) of lattice points in 
this hiperbox is given by the volume of the quotient between 
these lattices, i.e. 



Mi(t) 



det(Bt 



det(B) ' 

In this example Af 4 (1.35234) = 4.46213 x 10 32 . Since all 24 
torus layers are symmetric, the total number of point in the 
48-dimensional spherical code is given by 



M = 24 X Mi(t) = 1.07091 x 10 34 . 

Since the points inside the hiperbox V c are defined by a 
quotient of Abelian groups, we can use the standard Smith 
Normal Form [ 3 ) of the matrix Q = B~ 1 B\ to classify this 
group and to get the set of generators. 

In this example, in each layer the commutative group is 
isomorphic to Zg4 x Z256 x Z512 9 x Zn264- 

In addition, every point Xi G XXSC(48, 0.1), 
1 < i < 4.46213 x 10 32 in each layer of the spherical 
code can be generated as the orbit though a product of power 
of rotation matrices by a initial vector c(t) as follows 



x t = W{ 



rfei 



w* 2 



where 



.x , 



47 td>48 



x = (c(t)i,0,c(i) 2 ,0,--- ,c(*) 2 4,0) G S*' c 






< 




< 


63 





< 


k 2 


< 


255 





< 


kj 


< 


511 





< 


ho 


< 


11263 


1 


< 


3 


< 


9 



and Wi,W2,Wj,Wio represent, respectively, the 
generators of subgroups of rotation matrix in 
0(48) (orthogonal 48 x 48 matrices) isomorphic to 

2^64, Z256, Z512 9 , Zn264- 

Since the group structure are the same for all 24 layers in 
this TLSC(4&, 0.1) spherical code, there is a natural labeling 
for all the 1.07091 x 10 34 points in this code induced by the set 
of permutation vectors SC(L,d)+ and the commutative group 
code in each layer. It means that we are able to generated each 
one of the code points independently, which is a very useful 
property in many applications, specially channel coding and 
vector quantization. 

VI. Bounds and density of TLSC 

1 ) The grid TLSC: In this section we derive a lower and 
an upper bound for the number of points in a torus layer 
spherical code. Both bounds depend on a code SC(L, d) 
in L-dimensions. More specifically, to present the bounds 
we will assume that we have completed step (i) of the 
construction, i.e., assume we have selected and stored k points 
in SC{L,d)+. 

For given d, we construct a TLSC, by choosing Yx e as a 
subset of the rectangular lattice which lies in the hyperbox 

V c (pi. Let d = (cji )Ci2) ... )Ci£ ) G SC(L,d) + and u be a 
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point in the rectangular lattice is given by u = 2~2j=i tnjO>ij£ji 
where a,ij is the increment along the jth coordinate, 6j is the 
jth unit canonical basis vector and rrij is an integer, j = 
1,2, ...,L. Since the hyperbox V Ci has length 2wCij along 



Fig. 4. Illustration of a grid in a torus. 

the jth coordinate, we can determine the maximum number 
of lattice points in the hyperbox V Ci by finding the smallest 
an £ M. which satisfies 



d($ c (a tf ej), $c(0)) = 2cijsin ( , v 



> d. 



This leads to the following lower bound on the number of 
points in the code. 

Proposition 3: For a given SC(L, d) + with k points, we 
can design a TLSC(2L, d) with M(2L, d) points where 

k L 

M(2L,d)=J2l[W ij 

»=i j=i 

and 



arcsm 



2C,; 



if 



if 



2c, 



2a 



^ 1 



> 1 



2) An upper bound: Suppose again that we have k tori de- 
fined by a pre-selected SC(L, d) + . Let c, = (cq, c i2 , ■ • ■ , c^*) 
be the ith radius vector in SC(L,d) + . Without loss of ge- 
nerality we consider > cy+i Vi, j. We can use the bound 
presented in [ 1 1 1 [Proposition 7] to obtain the following upper 
bound for the number of points A/t c in each torus T Ci : 



Mr 



< 



(arcsin 3 ) L j=1 



L 

n r " A 



a code in that dimension. This process of projection should 
be carried out until a non zero value for is found. This 
is equivalent to place the points in a face of the hiperbox V Ci - 
We remark that, in the worst case all but one cy is zero and 
the torus degenerates to a circle where at least 



max 



2c it 



1 



points can be placed. 

Let Mj. , for some 1 < p < L, be the number of point that 
fit in the p-face of the hiperbox V Ci , (p = 1 corresponding to 
the one-dimensional degenerated torus obtained when just the 
first coordinate Cn is non zero an p = L corresponding to the 
i-dimensional flat torus obtained when all are non zero. 
Thus 



(arcsinf)^ 



il(cii)A ? 



Therefore the maximum number of points in each tori is 
given by 



M 



max 

i<p<fe 



Ml 



This allow us to dirive an upper bound for the number of 
points M(2L, d) of a TLSC(2L, d). 

Proposition 4: Given a SC(L,d)+ with k points, the num- 
ber of points in a TLSC(2L, d) satisfies 

k 

M(2k,d) < J2 M h- 
i=i 

Table |V] shows a comparison between these bounds and a 
TLSC(4, d) designed in the SectionllV] Note the tightness of 
the upper bound when the distance decreases. 



d 


TLSC(4,d) 


grid lower bound 


upper bound 


0.5 


172 


120 


194 


0.4 


308 


208 


360 


0.3 


798 


612 


826 


0.2 


2718 


2148 


2854 


0.1 


22,406 


18,884 


22,478 


0.01 


2.279 xl0 Y 


1.967 xlO 7 


2.279 xl0 Y 



where is the maximum center density of a packing in M L . 
However, there is an additional consideration. In case of the spherical codes 
smallest coordinate of Ci is small enough we may get 



TABLE V 

Bounds for 4-dimensional torus layer spherical codes at 
various minimum distances 



A. Density of TLSC 

In this section we analyze the density of the torus layer 



In this case, we should remove the smallest coordinate of Cj, 
to project the torus T Cj into one less dimension and consider 



Let T be the standard Gamma Function 

T(x)-- 

We shall use 



Sr 



Ltt l / 2 



r((L/2) + i) 



7 



for the (L— 1) -dimensional volume (surface area) of the sphere 

S^ 1 c R L and 



V L := 



r((i/2) + i) 



for the L-dimensional volume of the ball bounded by S L 1 . 
We also use 

SC(6/2,L) :=<Si_i / sirr L - 2 zcfa; 
Jo 

for the (L — 1) -dimensional volume of a spherical cap on the 
sphere S 21 "^ 1 with angular radius 8 = 2 arcsin(d/2). By using 
the Taylor series of sin(a;) and some standard calculations (see 
e.g. E) we can obtain 



SC(JB/2,L) = V L _! 



L-l 



0{d 



L+ls 



The density of a L-dimensional spherical code with mini- 
mum distance d and M codewords is proportion of the area 
of S L ~ 1 occupied by the union of the spherical caps centered 
at the codewords and with angular radius 9 = 2 axcsin(rf/2), 
that is, 

SC(d/2,L) M 



Asc = 



S L 



For a given minimum distance d, the maximum cardinality 
of a L-dimensional spherical code is unknown for all L > 3, 
except for a handful of values of d |[T4l . therefore the prob- 
lem of determine the maximum density of a L-dimensional 
spherical code is approached through bounds. 

Next proposition approaches the density of a torus layer 
spherical code, for asymptotically small d. 

In what follow we denote f(d) ~ g(d) if 



lim 



f(d) 



1. 



d->o g(d) 

Proposition 5: The torus layer spherical code density 
Atlsc is upper bounded and asymptotically approach the 
density of AA iX A L _ii where A„ is the densest lattice in R n . 
Proof: 

The torus layer spherical code density is given by 

SC(d/2,2L) M 



A 



TLSC 



S2L 



(7) 



where M = M, is the total number of codewords, Mj is 

i=l 

the number of codewords in the i-th torus layer and k is the 
total number of layers on which the code lays on. 

When the distance become small, the number of points 
in each layer can be approached by considering the best L- 
dimensional lattice packing fTD 



and therefore 



Mi 



M 



{d/2)^V L ' 



k 

E — 

i=l 



(d/2) L V L 



M 



k L 

L TY[2nc i 



(d/2)tV L ^l\ 



J2i-1 



can be foliated by flat tori, we may 



Since the sphere S 
assert 

k L 

i=l j=l 

where the element of (L — l)-volume dV is the volume of the 
"positive" part of the sphere S^ -1 divided by the number of 
tori, 



dV 



2 L 



^2 L SC(d/2,L-l) 

Therefore we may assert 



VTT(27r Cii ) 



i=i j=i 



SC(d/2, L - 1) 



S2L 
~dV'' 



and the number of codewords can be estimated by 

Aaz, S 2 l Aa z ,A Ai/ _ 1 5 2 l 



M 



(d/2) L V L sc(d/2,L-i) (d/2) L V L SC{d/2, L-l) 

A *L-1 



dct Ai, dct Al-i 



2LV- 



2L 



2LV- 



2L 



{d/2) L V L {d/2) L -W L ^ (d/2) 2i -idetA L detA L _ 1 
Thus, from Q, when d — > we get 

SC{d/2,2L)-. 



A 



TLSC 



' (d/2) 21 -- 1 dct A L dct Al - 



and so 



A 



S2L 
V2L-1 



TLSC 



det A L det A L _ 



which is the density of the Cartesian lattice A^ x A^_ x . ■ 
It should be remarked here that this asymptotic density is 
much better than the asymptotic density of apple peeling [6] 
construction but certainly worst than the best lattice packing 
density in i? 2L_1 , which can be achieved by the wrapped [8] 
and laminated [9] codes. On the other hand, for not that small 



d, as we have seen in Sec. IV a TLSC can outperform these 
previous constructions besides having the mentioned features 
inherit from its group structure. 

VII. Decoding 

Given an arbitrary x € R™ and a n-dimensional spherical 
code SC, the maximum-likelihood decoding problem is to find 



y = arg mm \\x - 

Vi€SC 

For any x £ E™ and a spherical code SC 



(8) 



arg mm \\x — y\\ 

y&SC 



arg mm 

yeSC 



IN 



In fact, let y g SC, such that 



x 

R 



< 



X 

Jx\ 



V zeSC. 
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Thus, 

2-2 
and 



y) <2-2 



z ) => (x,y) > (x,z) 



For most applications, we can conclude the decoding pro- 
cess assuming a suboptimal solution. We can also apply an 
additional step to get a maximum-likelihood decoding as 
follows. 

Let W£ G M L be the closest point to in V Ci and 



\\ x ~y\\ = IM| 2 +l-2 (x,y) < ||a;|| 2 +l-2(x,z) = ||ac-z||. Vi = il) Cf (w £ ) be its image in S 



2L-1 



If dp 



x < 



2' 



Therefore, for decoding any received vector x, using a 
spherical code, we can assume a; is a unit vector. 

One approach to solve <JSj is computing the all the inner 
products between x and y\ g SC and search for 



y 



arg max {x, y l 



This process requires 0(Mn + M) flops, but the main 
problem here is that this approach requires the storage of all 
codebook in the decoder, which is a restrictive requirement for 
many applications with limited memory. In addition, mostly 
computations required in the decoding process are done in the 
half of the code dimension. 

In what follows we address on decoding in a TLSC. 

For any 2L-dimensional unit vector x we can write 



7i 



X\ X2 

7i ' 7i 



x — I 7i cos — , sm — 
7i 7i 



,1l 



X2L-1 X2L 



8l . 0L 

. , 7i cos — , sm — 
JL JL 



Where, 



7i = 



1 < i < L 



X2i-1 

arccos I |7i 

7* ' 



1< i < L. 



This means that x belongs to a fiat torus of radius 
c x = (7i>72,-" ,7l)- In general, c x SC(L,d) + , i. e. c x 
does not defines a layer in the spherical code and we must 
project x in the closest layer. 

This process involves a spherical decoding in L-dimension, 
considering just the points in SC(L,d) + which defines the 
layers of tori in TLSC. As the number of tori is, in general, 
much smaller of code's cardinality, it does not increase the 
complexity of the entire process. 

For any Cj = (ca, c i2 , ■ • • , oil) € SC(L, d) + the vector 



Xi = | Cn \ cos — , sm — 
7i 7i 



. 9l ■ d A 

,c iL cos — ,sm — 

1L 1L J 



is the projection of x in the torus T Ci , i.e., 

\\x - x~i\\ < \\x - y\\Vy e T Ci . 

Let T C( be the closest torus to x. With high probability, the 
solution of (|8jl belongs to the torus T c? , and can be found by 
decoding the vector 



n = ^{xt) = y 



6>ic $1 6» 2 cj 2 



71 72 1L 

in the L-dimensional hyperbox V Ci using an efficient algo- 
rithm in the half of the code's dimension, which depends on 
the structure of the points in V cv For instance, for codes 
designed in previous section, the decoding in V Ci requires 
0(L) flops and does not need to store the codebook 



the maximum-likelihood decoding is over and j/j is the solution 
for fj. 

If d% > |, there might exist another w in some other torus 
T Ci such that 

\\ w ~ x \\ < \\Vi ~ x \\ 

Let us define precisely what tori must be checked. 
Let N = ^2) • • • the set of tori for which 



Ai = 



< d. 



6 



We will assume that A fj < A ?i+1 Vz = 1, 2, ...,j. 

Thus, we need to decode iteratively x in the torus defined 
in TV, getting a set of candidates Y = y^ , y& , ■ ■ ■ ,y^}, 
Y C TLSC(2k,d). 

The output of decoding will be the point y* E Y which 
satisfies 

\\y* -x\\<\\y- x\\ Vy e Y. 

In order to accelerate this process, each value of d^, obtained 
iteratively, can be used to reduce the set Af. Figure [5] illus- 
trates the decoding process in a TLSC(2L,d). Each circle 
represents a torus T c in the code. In this example, just the tori 
T c< and T Cw must be checked. 




Fig. 5. Decoding process in a TLSC(2L, d). 

The computational complexity in a 2L-dimensional spheri- 
cal codes constructed in layers of flat tori is dominated by the 
complexity of decoding in a L-dimensional hyperbox. 

VIII. Conclusion 

We propose a new construction of spherical codes based on 
the foliation of the unit sphere in even dimensions by fiat tori. 
Given a minimum distance d, the first step in this construction 
is to select torus layers which have minimum distance d. A 
codebook is then constructed in each layer by choosing a set 
of points in a hyperbox in half the code dimension. These 
points can be selected as cosets of a dense lattice in R L 



inducing a structured spherical code in R 8 which can be 
easily labeled and is generated by a commutative group of 
rotation matrix in each layer. The performance of these torus 
layer spherical codes is good when compared to the well- 
known wrapped spherical codes [8|, laminated spherical codes 
||9) and apple-peeling codes [6] for not asymptotically small 
distances. Concerning the coding and decoding process the 
main advantage comes from their homogeneous structure and 
the underlying lattice codebook in the half the code dimension. 
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